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Abstract 

The pseudoscalar-vector-vector (PVV) correlator is constructed using two meson multiplets in the vector and two in 
the pseudoscalar channel. The parameters are constrained by the operator product expansion at leading order where 
two or all three momenta are considered as large. Demanding in addition the Brodsky-Lepage limit one obtains (in 
the chiral limit) a pion-vector-vector (irVV) correlator with only one free parameter. The singly virtual pion transition 
form factor J r 7r o 77 » and the decay width of u> —> 7r°7 are independent of this parameter and can serve as cross-checks of 
the results. The free parameter is determined from a fit of the ui -f transition form factor . The resulting ttVV 

correlator is used to calculate the decay widths u> -A n°e + e~ and u> -A- n W - and finally the widths of the rare decay 
7r° —> e + e~ and of the Dalitz decay 7r° —> e + e _ 7- Incorporating radiative QED corrections the calculations of 7r° decays 
are compared to the KTeV results. We find a deviation of 2 a or less for the rare pion decay. 

Keywords: 13.20.Cz Decays of 7r mesons, 13.40.Gp Electromagnetic form factors, 11.30.Rd Chiral symmetry, 12.40.Vv 
Vector-meson dominance 


1. Introduction and summary 

Two major challenges of contemporary particle physics are 
the search for beyond-Standard Model physics and a bet¬ 
ter understanding of the non-perturbative low-energy sec¬ 
tor of the strong interaction. Typically both aspects inter¬ 
mix when it comes to high-precision determinations of low- 
energy quantities and the corresponding Standard Model 
predictions. If a low-energy observable is potentially in¬ 
fluenced by quantum effects from new particles, then it 
is also influenced by hadronic loop effects. The latter of¬ 
ten — if not always — constitute the main uncertainty 
of the Standard Model prediction. On the other hand, 
new physics can only be revealed if a significant devia¬ 
tion between experiment and Standard Model calculation 
is observed. Of course, this requires small uncertainties 
for both the experimental and the Standard Model result. 
Our poor understanding of non-perturbative QCD could 
provide a serious hurdle for a reliable Standard Model cal¬ 
culation and/or for a reliable uncertainty estimate of such 
a calculation. 

Two quantities of current interest which might indi¬ 
cate some deviation between experiment and the Stan¬ 
dard Model are the gyromagnetic ratio of the muon [1, 2] 
and the rare decay of the neutral pion into electron and 
positron [3, 4], An important quantity that enters both ob- 
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servables is the pion transition form factor, i.e. the three- 
point correlator between a neutral pion and two electro¬ 
magnetic currents. 

Two tasks are carried out in the present work: First, the 
pion-vector-vector (nVV) correlator is determined by com¬ 
bining high-energy quark-based information with low(er)- 
energy hadronic information. We follow the general ap¬ 
proach proposed in [5] with some refinements. Second, 
we explore some phenomenological consequences of our 
correlator and focus in particular on the rare pion decay 
7T° —>• e + e - . Here we include also QED radiative correc¬ 
tions along the lines of [6, 7]. 

Two limits of QCD are of particular interest for low- 
energy hadron physics: The chiral limit where the masses 
of the two or three lightest quarks are neglected [8, 9, 10, 
11, 12] and the limit where the number of quark colors, N c , 
is sent to infinity [13, 14]. Concerning the nW correlator 
the chiral anomaly fixes the low-energy strength unam¬ 
biguously in the chiral limit. For a large number of col¬ 
ors there are infinitely many, infinitely narrow, i.e. stable, 
quark-antiquark states for every combination of quantum 
numbers. They show up as poles in the n-point correlators 
of quark currents. In the real world of three colors the 
hadrons generically turn to unstable resonances because 
the hadronic interactions do not vanish any more. The 
poles in the correlators turn to cuts (and poles in other Rie- 
mann sheets). The cuts start at the corresponding rnany- 
body thresholds. Thus the relevance of the large- N c limit 
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for the real world is highest, if one considers hadrons which 
are narrow and/or kinematical regions where there are no 
(significant) cuts. This is our guiding principle when ex¬ 
ploring phenomenological consequences. Concerning form 
factors there are no cuts in the space-like region. 

In the Standard Model the rare decay of the pion into 
electron and positron is caused by a loop where the pion 
first turns into a pair of (real or virtual) photons; see 
Fig. 4 below. At this vertex the pion transition form fac¬ 
tor sneaks in. If the form factor was replaced by a con¬ 
stant, the loop would diverge [15, 16, 17]. In QCD the 
pion transition form factor is suppressed for large virtu¬ 
alities [18, 19, 20]. This leads to a finite result for the 
7r°-e + -e _ amplitude at the one-loop level of the Standard 
Model calculation. Thus for a quantitative determination 
of the branching ratio of the considered rare pion decay it 
is necessary to know where and how fast the pion transi¬ 
tion form factor reaches its asymptotic form which in turn 
depends on the various combinations of virtualities. 

These considerations show that one needs information 
from various QCD regimes: The threshold regime gov¬ 
erned by the chiral anomaly, the regime of hadronic reso¬ 
nances, and, finally, the regime of asymptotically high en¬ 
ergies dictated by quarks and asymptotic freedom. These 
regimes are connected in the approach of [5] where the op¬ 
erator product expansion (OPE) for various three-point 
correlators of quark currents is matched to a hadronic 
ansatz that satisfies the chiral constraints for the low- 
energy limit. This ansatz is furthermore based on a trun¬ 
cation of the infinite tower of stable hadronic states that 
appears in the limit of a large number of colors. 

In principle one can work out arbitrary many orders in 
the OPE and match to the parameters that emerge with 
the tower of hadron states. However, the higher orders in 
the OPE contain unknown quark and gluon condensates 
of high dimensionality. Thus in practice the model de¬ 
pendence emerges from a selection of the to be matched 
OPE constraints and from the choice where to truncate 
the tower of hadron states. Using one hadron multiplet 
per channel and leading-order OPE constraints has been 
studied in detail in [5]. This truncation is called “lowest- 
meson dominance” (LMD). In the present work we will 
explore the consequences of having two hadron multiplets 
per channel. In the language of [5] our approach would 
be called “LMD+V+P”. To avoid this clumsy name we 
decided to introduce the name “two-hadron saturation” 
(THS). 

Concerning our quantity of interest, the ttVV corre¬ 
lator and the corresponding pion transition form factor, 
the starting point on the level of quark currents is the 
pseudoscalar-vector-vector ( PVV) correlator. The con¬ 
sequences of LMD for this quantity have been studied 
in [5]. An application to the rare pion decay to electron 
and positron was presented in [17]. Two vector multiplets 
have also been considered in [5, 21]; see also [4] where this 
has been used for the rare pion decay. What makes our 
approach different from previous works is that we explore 


in detail the consequence of two multiplets in each channel 
and that we fit and/or compare to data on the 7r°wU cor¬ 
relator. In fact, including a second multiplet in the vector 
channel involves the energy region of about 1.4 GeV [22]. 
In this region there is also a pseudoscalar multiplet. Thus 
the extension from LMD to two multiplets for a channel 
suggests to use two multiplets for every channel. Con¬ 
cerning the second aspect, the interrelation to the 7 t°loV 
correlator, we will come back to this issue below, after 
discussing in more detail the pertinent high-energy con¬ 
straints. 

There is yet one more short-distance limit to be consid¬ 
ered. Instead of studying the high-energy limit of correla¬ 
tors of quark currents (OPE), one can also study the high- 
energy limit of correlators that involve specific asymptotic 
states like hadrons or photons together with one or sev¬ 
eral quark currents. In particular the high-energy behavior 
of the pion-photon-vector correlator has recently gained 
much attention since the BaBar data [23] seem to contra¬ 
dict the Brodsky-Lepage (B-L) scaling limit [20] while the 
Belle data [24] seem to support it. 

Using THS we are able to satisfy all leading-order OPE 
constraints for the PVV correlator and in addition the B- 
L constraint for the pion-photon-vector correlator. While 
LMD satisfies the same OPE constraints, it violates the 
B-L constraint as can be easily deduced from the explicit 
form given in [17]; see also the discussion in [5]. As we 
will show below, the constraints from the leading order 
of the OPE, from B-L and from the chiral anomaly to¬ 
gether fix the THS approach to the nVV correlator up 
to one single parameter, which we call n. If the invariant 
mass of one of the vector currents in this correlator is set 
to zero, i.e. for the pion-photon-vector correlator, then n 
drops out. In other words we have full predictive power 
for this correlator. Aiming at the rare pion decay into 
electron and positron one needs the full information on 
the 7 tVV correlator for arbitrary invariant masses of the 
two vector currents. In principle, the parameter k would 
be best determined from data on the ttVV correlator with 
both invariant masses being different from real photons. 
Unfortunately such data do not exist. 

In this situation we turn to the second-best choice. Pro¬ 
jecting the 7 tVV correlator on one of the vector mesons 
that we include in THS yields a 3-point correlator for the 
pion, the vector meson and a quark current with vector 
quantum numbers. Given that our approach is based on 
the large-A^ c limit where mesons are approximated by in¬ 
finitely narrow states, it is suggestive to use a narrow vec¬ 
tor meson. Since phenomenologically and in the large- N c 
limit the pion decouples from the meson [25, 26, 27, 14], 
we are left with the w meson as the best choice for a vector 
meson. Consequently we will use data on the uj-n transi¬ 
tion form factor [28] to fix our remaining parameter n. 

The 7rUU correlator obtained in this way from THS 
shows an awesome behavior: If the virtuality of one vec¬ 
tor current becomes large while the other is set to zero 
(photon case), the asymptotic B-L limit is reached rather 
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fast, resembling essentially strict vector-meson dominance 
(VMD) [29, 30]. The scale that defines where the approach 
to the asymptotic limit sets in is basically given by the 
mass of the p/ui meson. On the other hand, if both vector 
currents have the same large virtuality, the corresponding 
asymptotic limit is reached very late for the ttVV correla¬ 
tor as obtained from THS. This finding points to the rele¬ 
vance of details of hadronic physics above 1 GeV. For this 
case the ttVV correlator from VMD falls off much faster 
than demanded by QCD, while for LMD the asymptotic 
limit is reached much earlier than for THS. Since the rare 
decay n 0 —>■ e + e~ is sensitive to both high- and low-energy 
physics, it is interesting to study how this intriguing be¬ 
havior of the ttVV correlator obtained from THS influ¬ 
ences the branching ratio of this rare process. 

Before going into the details of our findings we shall 
compare our approach to related ones from the literature. 
The uj-tt transition form factor and related quantities have 
also been addressed in [31, 32, 33] where the ground-state 
vector mesons are treated as light degrees of freedom. By 
construction the approach is restricted to low energies, i.e. 
high-energy constraints were not considered. Nonetheless, 
it turns out that the THS result for the uj-tt transition 
form factor is numerically very close to the one from [31]. 

Conceptually close in spirit to THS is the Lagrangian 
approach utilized, e.g., in [34, 35] (earlier references can 
be traced back from these works). Also here hadron res¬ 
onances in the large- N c limit are used to interpolate be¬ 
tween the low-energy region governed by chiral perturba¬ 
tion theory and the high-energy region governed by the 
OPE or quark scaling considerations. In [35] a PVV and 
a corresponding ttVV correlator are constructed with one 
multiplet in the vector channel and two multiplcts in the 
pseudoscalar channel (i.e. the Goldstone bosons and one 
resonance multiplet). For an extension to two multiplets in 
the vector channel see [36] . The ttVV correlator from [35] 
satisfies the B-L constraint, but the PVV correlator satis¬ 
fies only one leading-order OPE constraint, not all of them. 
As shown in [37] the Lagrangian utilized in [35] is not capa¬ 
ble to provide correlators that satisfy all OPE constraints; 
see also the discussion in [5] . Instead of the other OPE con¬ 
straints on the PVV correlator a high-energy constraint 
based on quark counting rules is imposed on the npV cor¬ 
relator in [35]. 

In our work we impose high-energy constraints on the 
PVV correlator (THS satisfies all leading-order OPE con¬ 
straints) and on the nVV correlator (the B-L limit), but 
not on the tt pV or 7 tujV correlator. Our philosophy is that 
we take the first multiplets to resemble the correspond¬ 
ing ground-state physical particles but the second multi¬ 
plets to mimic the effect of the tower of infinitely many 
excited states. Since we study with PVV an order pa¬ 
rameter of chiral symmetry breaking [5] one expects that 
the second multiplets are close to the physical states that 
are the first excitations on top of the ground states. How¬ 
ever, the weighted average of the whole tower of states 
might shift the effective mass higher up. In the present 


work we explore the uncertainty of the THS approxima¬ 
tion by changing the masses of the second multiplets from 
the first to the second physical excitations. Concerning 
the 7 tVV correlator we expect to obtain reasonable results 
because the pion might not resolve too many details of the 
intermediate-energy region. Thus replacing the tower of 
excited resonances by the respective lowest excitation and 
demanding high-energy constraints for the ttVV correla¬ 
tor might be good enough. In contrast, a 7r pV or ttujV 
correlator resolves more from the intermediate-energy re¬ 
gion because the vector meson induces a larger mass scale. 
(The same would apply to a pPV correlator.) Therefore 
to satisfy high-energy constraints in this case we expect 
that one would need a more detailed modeling than just 
having THS, i.e. one excitation (plus the ground state) 
in each channel. It might be worth to explore a three- 
hadron saturation scenario, but this is beyond the scope 
of the present work. Therefore we demand constraints on 
the PVV and ttVV correlators, but disregard constraints 
for three-point correlators of one pseudoscalar meson, one 
vector meson and one vector quark current. For the same 
reason we only consider leading, but not subleading high- 
energy/OPE constraints. 

From a formal point of view our approach is close to the 
successive Pade approximations as utilized, e.g., in [38]. 
Our correlators are also approximated by rational func¬ 
tions. In our approach, however, the poles of the correla¬ 
tors are related to physical states (in the large- N c limit). 
In contrast, in the Pade framework of [38] one determines 
the rational functions by fits to data. It is not the pur¬ 
pose of this Pade approach to look for the poles of the 
obtained rational functions. In fact, there are no physical 
restrictions from outside that would make sure that these 
poles correspond in any way to physical hadrons. But this 
is what S-matrix theory suggests (in the large- N c limit): 
Anything beyond polynomials should be caused by unitar- 
ity, analyticity, crossing symmetry and physical states. 

Our work is complementary to the dispersion theoretical 
approach of [39, 40]; see also [41]. While we cannot expect 
to reach the accuracy of a dispersive approach concern¬ 
ing low-energy quantities, our framework has the advan¬ 
tage to provide a smooth and physical connection between 
the low- and high-energy region and between the quark 
based and hadron based correlators. In practice, disper¬ 
sive calculations are based on an excellent account of the 
low-energy region (up to about 1 GeV) and a high-energy 
completion, i.e. a matching to the high-energy behavior 
deduced from S-matrix theory, QCD or QCD related ap¬ 
proaches; see the discussion in [40] concerning the ttVV 
correlator. As already discussed, for our doubly virtual 
nVV correlator there are regions in the virtualities where 
the asymptotic regime is reached only at very high ener¬ 
gies. In practice this might imply that a naive matching 
of a dispersive calculation to the asymptotic regime might 
miss part of the physics present at intermediate energies. 
Clearly it is worth to explore the interplay of a dispersive 
calculation with THS in the future. 
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The rest of the paper is structured in the following way: 
In Section 2 we construct the PVV and nVV correlators 
subject to high- and low-energy constraints. The results 
for the nVV correlator are compared to data and to other 
approaches (LMD and VMD). It is studied how the shape 
of the correlator changes when varying the remaining free 
parameter k in a reasonable range. In addition, a model 
uncertainty is estimated by varying the mass of the second 
vector-meson multiplet between the physical masses of the 
first and second excitation. 

In Section 3 the iuxV correlator is constructed. The 
parameter k is determined from a fit to the lu-tt transi¬ 
tion form factor. The widths or branching ratios, respec¬ 
tively, for the corresponding decays w —» 7r°7,7r°e + e _ , and 
7 r°/i + p,“ are determined for THS (and also for LMD and 
VMD). Actually the first decay does not depend on k. We 
find good agreement between THS and the experimental 
results. 

We address the rare pion decay to electron and positron 
in Section 4. Including radiative corrections the branching 
ratio of this process is calculated and we compare again 
the THS result to other approaches. For the THS case 
this branching ratio is sensitive to k. Direct comparison 
to the experimental value from KTeV [3] seems to suggest 
that a discrepancy persists on the level of 2 a. 

In Section 5 we study the properties of the singly vir¬ 
tual pion transition form factor in the low-energy region, 
such as slope and curvature. Note that these quantities 
do not depend on the parameter n. For the Dalitz decay 
7T° —>• e + e _ 7 we calculate the decay width taking into ac¬ 
count next-to-leading order (NLO) radiative corrections. 
These are evaluated along the lines of [42] . Again we com¬ 
pare THS to other approaches. With the full set of radia¬ 
tive corrections at hand, we take a fresh look on the KTeV 
result. Considering some radiative corrections that were 
not accounted for in the analysis suggests that the previ¬ 
ously stated discrepancy might be even reduced to 1.5 a. 
The main message here, however, is that the radiative cor¬ 
rections are now theoretically under control [6, 7, 42] and 
can be used in future data analyses. 

In Section 6 we provide an outlook how THS can be 
further utilized and how the scheme can be extended. Fi¬ 
nally, an appendix is added to provide the pseudoscalar 
form factors in terms of the loop integrals required for the 
rare pion decay. 

Comparing in detail various approaches throughout the 
present work reveals that the VMD form factor proves 
to be phenomenologically very successful in most applica¬ 
tions, in spite of the facts that VMD is so simple and par¬ 
tially possesses an improper high-energy behavior. How¬ 
ever, the THS form factor, which satisfies all the consid¬ 
ered constraints, works very well, too. Still, the mean 
value of the KTeV result remains a challenge for all ap¬ 
proaches, even though the discrepancy has been reduced 
by a considerable level. We provide a brief discussion on 
this point in the corresponding sections. 


2. THS approach to the 7 rVV correlator 

Following [5, 37] we introduce the PVV correlator 
H (r 2 * ;p 2 ,q 2 ) by 

d abc e liVol pp o ‘q 0 Ii.{r 1 ] p 2 , q 2 ) 

r (1) 

= / d A xd 4 ye ip ' x+tq ' y (0\T[P a (0)V b {x)V v c (y)}\0) 

with r = p + q. The vector and pseudoscalar current, 
respectively, are defined by 1 


V“{x) = q(x)'y ll T a q(x) , P a { x) = q(x)ij 5 T a q(x ). (2) 

In the above formulae we have used 

Tr [T a ,T b ] = i<5 ab , d abc = 2Tr [{T a ,T b }T c ] . (3) 

For a = 1,..., 8 we have T a = X a /2 where A a denote the 
Gell-Mann matrices in flavor space. Since we utilize the 
large- N c limit we have to deal with flavor nonets, i.e. with 
U(3) instead of SU(3). The formulae (3) provide a natural 
extension to a = 0. 

Working in the chiral limit 2 and considering two me¬ 
son multiplets in the vector and two in the pseudoscalar 
channel, we propose a correlator of the form 


n THS (r 2 ;p 2 , q 2 ) = 

X {p 2 - My ± ) (p 2 


B 0 F 2 

r 2(r 2 _ m 2 ) 

P(r 2 -p 2 ,q 2 ) 

‘ My 2 )(q 2 — My i ){q 2 



( 4 ) 


Here P(r 2 ;p 2 , q 2 ) is the most general polynomial in its ar¬ 
guments that is symmetric in the second and the third 
argument. For our purpose, a polynomial of order four 
is sufficient inasmuch as higher powers of momenta in the 
numerator are not allowed due to the desired high-energy 
behavior that we will specify below. Thus we start with a 
term containing 22 monomials and associated free param¬ 
eters. Schematically this looks as follows: 

P(r 2 ;p 2 , g 2 ) = c 0 pV + d[(p 2 ) 4 + (g 2 ) 4 ] 

+ c 2 [(p 2 ) V + (T)V] + c 3 (r 2 )Vg 2 + ... . ( } 


The quantities F and B 0 in (4) denote the usual low- 
energy constants of chiral perturbation theory [11, 12]. We 
will specify them further when needed. Mp denotes the 
mass of the second multiplet, i.e. first excitation, in the 
pseudoscalar channel. The first multiplet (ground state) 
is, of course, the massless multiplet of Goldstone bosons. 

The masses of the lowest two vector-meson multiplets 
are denoted by My 1 and My 2 . As already spelled out in 


1 Our convention is 75 = vVA 

2 The chiral limit is only used for the construction of the correlator. 

Once the form factor is settled, we take the physical pion mass for all 
the kinematics used for comparing the predictions to experimental 

data. 
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Section 1 we use for My 1 the mass of the ground-state 
vector-meson multiplet (in the chiral limit). We assume 
that this mass is approximately given by the mass of the 
oj or p meson. For My 2 it is suggestive to use the physical 
mass of the first excitation in the vector channel. However, 
we do not use a fixed mass here, but study the impact of a 
variation of My 2 on our results. In this way we explore the 
uncertainty caused by the higher-lying excitations that are 
neglected in THS. In practice we vary My 2 in the range 

My 2 e [1400,1740] MeV, (6) 


which is the interval between the masses of the first and 
the second physical excitation [22]. 

Finally we note that the same logic applies to the mass 
Mp of the second pseudoscalar multiplet. For the PVV 
correlator one should study the impact of a variation of 
Mp. As we will see below, however, this mass does not 
show up in the final expression for the tyVV correlator. 

We now demand that the ansatz (4) satisfies all the rel¬ 
evant high- as well as low-energy constraints in order to 
minimize the number of free parameters introduced in (5). 
Starting with the general correlator (4) we apply the fol¬ 
lowing leading-order OPE constraints [5]: 


n((Ar) 2 ; (Ap) 2 , (Xq) 2 


2 B ° F A 4 


1 r 2 + p 2 + q 2 


r 2 p 2 q 2 


O 


1 


n(r 2 ; (Ap) 2 , (r - Ap) 2 ) = B 0 F 2 -+ O 
It turns out that the third OPE constraint [5] 


n((<? + Ap) 2 ; (Ap) 2 ,g 2 ) = ^^/(9 2 ) + e>(^ , 


(7) 

( 8 ) 

(9) 


is automatically fulfilled by our ansatz. Here / denotes a 
function (actually a two-point correlator [5]) that depends 
only on q 2 and not on the other kinematic variables p 2 and 
Pi- 

Next, we define the ttVV correlator: 


Fnvv{p 2 ,q 2 ) = lim 7 

Ar r 2 —>0 


where 

= ^(0|(u7 5 u - J75gO|7t 0 ) (11) 

denotes the overlap between the pion field and the pseu¬ 
doscalar quark current. With the usual conventions from 
chiral perturbation theory one obtains Z = BqF [11, 12]. 
For large A one then finds 


•7wv((Ap) 2 , (Ap) 2 ) 

~ lim r 2 n(r 2 ; (Ap) 2 , (r - Ap) 2 ) 

r 2 —>0 

b 0 f b ° f a 2 p 2 


( 12 ) 


o 


A 3 


which means that in agreement with [18, 19, 5] we have 

F 

Fnwiq W) -t -o , q 2 -> -OO. (13) 

q 2 


Hence, this condition is satisfied automatically on ac¬ 
count of the OPE constraints. Actually, for the quantity 
Fnvviq 2 , q 2 ) also the subleading order in the high-energy 
expansion is known in terms of a quark-gluon condensate 
[43]; see also the corresponding discussion in [1]. How¬ 
ever, we refrain from incorporating this as an additional 
constraint. The reason is the same as to why we do not 
use constraints from transition form factors of hadronic 
resonance states. One might become too sensitive to the 
details of the intermediate mass region where we use one 
hadronic state to describe effectively the whole infinite 
tower of large-JV c excited states. If one used three in¬ 
stead of two hadronic states per channel, the subleading 
high-energy constraints might provide a viable input to pin 
down the growing number of resonance parameters. But 
this is clearly beyond the scope of the present work. 

Instead of involving subleading orders in the high-energy 
expansion we apply the B-L constraint [20] 


TnVvjO ,q 2 ) 24tt 2 F 2 1 . 

•7W(0,0) N c q 2 ’ q 


(14) 


We define the pion transition form factor as 

Tv° 7 * 7 *(pW) = ^7rvv{p 2 ,q 2 ) (15) 

and match at the photon point to the chiral anomaly, i.e. 
to the Wess-Zumino-Witten term [9, 10], 

3 N 

-^w(0,0) = -j;o 7 , 7 ,(0,0) = -^. (16) 

Together the constraints (7), (8), (14) and (16) provide 
us with a ttVV correlator that appears to have only one 
free dimensionless parameter n: 


J-™s ( p W) = _ 


N c 
8 t t 2 F 


m^m 4 2 


(P 2 - My i )(p 2 - My 2 ) {q 2 - My i )(q 2 ~ M £ ) 
4 :ir 2 F 2 (p 2 + q 2 ) 


1 + 


K p 2 q 2 


2 N c (4t tF) a 


n c m 2 m 2 2 


6 + 


2 2 
p q 


m$m 2 2 , 


(17) 


Note that our result is independent of the mass Mp of the 
first pseudoscalar excitation. This happens due to the fact 
that at the end of the day we could conveniently rescale 
the only free parameter left. From the structure of the 
result (17) it can be read off that k emerges from cq in (5) 
on account of 


co 


KMpMy^My 2 

(4t xFf 


(18) 


We note in passing that k scales with IV 3 . 
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A comparison to the work of [5] is in order here concern¬ 
ing the number of free parameters. We shall compare our 
THS approach to LMD+V of [5], i.e. to the case of one 
pseudoscalar and two vector multiplets. We recall that 
THS = LMD+V+P in the language of [5]. After applying 
the OPE constraints to THS we are left with 12 free pa¬ 
rameters for the PVV correlator. This compares to 7 free 
parameters for the case of LMD+V. The low-energy con¬ 
straint (16) fixes always one more parameter. Once we fo¬ 
cus on the 7r VV correlator we are left with 4 parameters for 
THS. For LMD+V one has 3. Demanding that the 7 rVV 
correlator of (14) drops like 1/q 2 limits the free param¬ 
eters to 2 for both cases THS and LMD+V. Demanding 
(14) quantitatively and not just a scaling with 1/q 2 yields 
one free parameter for THS as well as for LMD+V. To 
summarize, after applying all constraints THS has more 
parameters for the PVV correlator than LMD+V: 8 in 
THS versus 4 in LMD+V. But concerning the nVV corre¬ 
lator one ends up with one free parameter in both cases. 

According to [4] the 7 rVV correlator should satisfy the 
inequality 

\^vv{q\q 2 )\ < \3~ irvv(0, q 2 )\, q 2 < 0 . (19) 


It turns out that the THS expression (17) satisfies (19) for 
—45 < k < 30. We will see below that the values for n 
that we obtain from fitting to experimental data he well 
within this range. 

Before further constraining n from data it is illuminating 
to study the qualitative shape of the 7 rVV correlator when 
k is varied. In addition, we compare our results to sim¬ 
ilar approaches from the literature, namely to the VMD 
correlator [29, 30, 5, 37] 


^(pW) 


N c _ 

8n 2 F (jp 2 - M2.)(g2 - M 2 J 


( 20 ) 


and to the LMD expression [17, 5] 


x-LMD 

-rVV 


(pW) = j-^V>W) 



47 T 2 F 2 (p 2 + q 2 )' 

^ . ' 
( 21 ) 


See also [44] for a comparison of various correlators. 

First of all, we note that the VMD result grossly violates 

(13) while it is not far off from the B-L constraint (14). On 
the other hand, LMD satisfies (13), but does not satisfy 

(14) . By construction the THS correlator satisfies all the 
mentioned constraints, but it is interesting to see how fast 
or slow the asymptotic limits are reached. This is shown in 
Fig. 1. To facilitate a comparison with data we show the 
pion transition form factor (times the virtuality q 2 ) instead 
of the 7rW correlator. Note that the relation (15) between 
the pion transition form factor and the ttVV correlator 
only amounts to a rescaling. 

In the first panel of Fig. 1 we display the symmetric 
doubly virtual pion transition form factor. This plot shows 
three different types of lines and a gray band. We shall first 




- q 2 [GeV 2 ] 


Figure 1: Symmetric doubly virtual (first panel) and singly virtual 
(second panel) pion transition form factor as a function of the vir¬ 
tuality in the space-like region. The gray bands constitute the final 
THS predictions. See the main text for details how the uncertainties 
are determined. In the second panel the gray band nearly collapses 
to one full line. The thin full lines in the first panel show what hap¬ 
pens if k, is varied in a large range. The corresponding curve that 
constitutes the center of the gray band has been made bold. The re¬ 
sults from VMD and LMD are given by the dashed and dash-dotted 
lines, respectively. Data are taken from [45, 23, 24], 


discuss the lines: The dash-dotted LMD line approaches 
very quickly the asymptotic QCD result given by (13). 
The dashed VMD line falls stronger than what is required 
by QCD. The full lines labeled by values for n show the 
THS result for a mass of My 2 taken in the middle of the 
interval (6). All the full lines will approach the LMD line 
at very large momenta. One sees, however, that typically 
the THS lines reach this limit rather late. In particular, 
for k > 15 the THS form factor peaks at rather high mo¬ 
mentum values and with large magnitude in comparison 
to VMD and LMD. It reaches the asymptotic limit only 
very slowly. For the case K~9we see that the associated 
line is very close to LMD, but peaks later and approaches 
the asymptotic limit from above. For k ~ 6 we get the 
limiting case where THS is monotonically growing. Going 
further down with k a wiggle appears, which for k = 0 and 
low virtualities resembles very much the VMD behavior. 
For negative n the THS results start to undershoot VMD 
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at low virtualities. Moving to even lower re values there is 
always a region where the THS result becomes negative. 
Intuitively we find it hard to believe that the transition 
form factor of a ground-state hadron (the pion) would dis¬ 
play so many wiggles as suggested by the full lines for 
small values of re. Indeed, the determination of re, which 
will be carried out in the next section, reveals a re value 
of about 21. Thus we find the qualitative situation that 
the THS result for the pion transition form factor (times 
the virtuality q 2 ) overshoots the asymptotic limit, peaks 
at rather large momenta and approaches the asymptotic 
limit rather slowly from above. 

The gray bands in both panels of Fig. 1 constitute our 
final THS predictions for the pion transition form factor. 
The results from the next section are anticipated where re 
is further constrained. To obtain the gray bands all input 
for our THS pion transition form factor (F, My 1 ) is varied 
within the respective experimentally allowed regions, re = 
21 ± 3 as described in the next section, and My 2 is varied 
within region (6). The whole gray band in the first panel 
shows the qualitative behavior described previously. It 
also points to a significant quantitative uncertainty of our 
prediction for the symmetric doubly virtual pion transition 
form factor. In other words, data on this form factor or in 
general on any doubly virtual pion transition form factor 
would be highly welcome to better constrain our approach. 

In the second panel of Fig. 1 we display the singly virtual 
pion transition form factor. This plot shows two lines, a 
very narrow gray band and data. We first discuss the gray 
band, which nearly resembles a full line since its width is so 
small. Obviously the narrow gray band of the second panel 
contrasts with the rather broad band of the first panel. 
However, the singly virtual pion transition form factor is 
independent of re as can be easily deduced from (17) when 
putting p 2 to zero. The largest uncertainty comes from a 
variation of My 2 according to (6). Even this variation does 
not lead to very different results, which causes the curves 
to nearly collapse to one single line. Thus THS has rather 
high predictive (or rather postdictive) power for the singly 
virtual pion transition form factor. Another qualitative 
difference between the gray bands of the first and the sec¬ 
ond panel is the fact that the THS result for the singly 
virtual pion transition form factor reaches the asymptotic 
(B-L) limit rather early. We regard this intriguing behav¬ 
ior of our correlator as one of the highlights of our work: 
The THS pion transition form factor shows an early (late) 
onset of the asymptotic behavior for the singly (symmetric 
doubly) virtual pion transition form factor. It would be 
extremely interesting to see if this is supported by data in 
the future. 

Finally, we compare the THS result for the singly vir¬ 
tual pion transition form factor to VMD and LMD and 
to data. The LMD result for J r n .o 7 . 7 . (0, q 2 ) behaves as a 
constant for large \q 2 \. Thus, it cannot even qualitatively 
explain the data since, of course, it diverges after being 
multiplied by q 2 . This behavior can be seen in the sec¬ 
ond panel of Fig. 1 where the LMD result is given by the 


dash-dotted line. The results for VMD (dashed line) and 
THS are fairly close. CLEO data [45] are only available 
for low virtualities (— q 2 < 8GeV 2 ) and prefer VMD to 
some extent though THS appears to be acceptable as well. 
The situation turns around for higher virtualities. THS 
describes the Belle data [24] better than VMD. The early 
onset of the asymptotic behavior of the THS result for the 
singly virtual pion transition form factor makes THS es¬ 
sentially identical to the B-L limit. Not surprisingly then, 
THS agrees well with the Belle data and is at odds with 
the BaBar data [23] that are not fully compatible with the 
Belle data. We do not attempt to contribute to a clarifica¬ 
tion of the differences between Belle and BaBar. In view 
of these complications we have decided to stick to the B-L 
constraint (14) right from the start. 


3. Phenomenology of u decays 

Our strategy now is to obtain the parameter re of the newly 
proposed ttVV correlator (17) from the ttujV correlator. 
Thus we turn to ui data for the reasons specified in Sec¬ 
tion 1, namely the lack of doubly virtual data for the pion 
transition form factor and the fact that the ui meson is 
fairly long lived to resemble the situation of the large- N c 
approximation. 

We introduce the overlap between an ui meson and the 
vector current V, i.e. 

A w ) = i(0|(u7^u + d'y l _ l d)\uj{p, Xu)) , (22) 

where we assume that the u> meson does not contain hidden 
strangeness, which is a fairly good approximation to the 
real world [22]. For later convenience we also define the 
7 -w coupling strength 


Fui 


Zu 

Mu ’ 


(23) 


the modulus of which agrees with Fy as introduced in [46] . 

With this at hand we can obtain Fu from the u> —>• 
e^e - decay process. A direct calculation from the Lorentz 
invariant matrix element 


iM u 


ieZu 


.(-i) 


e+e~ = —jj— w(9l)(-*e7^)u(g 2 )-^2-c M 


(p, Au), 


(24) 


yields after averaging and summing over polarizations 


F u , = 


47T« 2 


27Mu F(w -> e+e-) 


j8e 


2 TO 2 

~M* 


(25) 


where p e u = \J 1 — Am1/M% is the speed of the electron 
in the rest frame of the decaying u> meson. Taking into 
account that m 2 <C M 2 , we can write 

Fu ~ \108n MuT (lo ^ e+e~). (26) 

e z 
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Using the values B(lo —> e + e ) = (7.28 ±0.14) x 10 5 and 
T(w) = (8.49 ± 0.08) MeV [22] we find T(w ->■ e+e") = 
(0.62±0.02) keV. Together with M u = (782.65±0.12) MeV 
we obtain from (26) the following value for the coupling 
strength: 

F u = (140 ± 2) MeV. (27) 

It is convenient to introduce the 7r°u;V transition form 
factor by 

•TVwO? 2 ) = -k~ hm (28) 

p 2 — 


This turns out to be the central quantity to determine k 
and to address the phenomenologically interesting decays 
oj —>• 7r°7 and oj —>• n°£ + £~ , where (. denotes a lepton. 

For the special case when q 2 = 0 we find (see also [34] for 
a similar expression for the case of one vector multiplet) 


x-ths 

•'V’wV 


( 0 ) 


1 N c M 2 M\ r 247r 2 F 2 " 
Z u (My 2 — My^j [ N c M 2 2 _ 


(29) 


Hence the decay width for w —> 7r°7 is independent of 
k, which provides us with full predictive power for this 
particular decay and can therefore serve as a cross-check 
of our formalism. We will come back to this decay after 
having determined k from the measured lo-tt transition 
form factor. 

In fact, the parameter k, which constitutes the only un¬ 
known parameter of our proposed 7 rVV correlator, can be 
obtained by a fit to the NA60 data [28] for the normalized 
transition form factor 




7 r°LdV 




*^ 7 T°Ujv(.Q ) 

•Ftt°ujV (0) 


(30) 


in the low-energy time-like region of q 2 with the result 


k = 21±3. 


(31) 


The result is displayed in Fig. 2 as a gray band. The 
uncertainty comes mainly from the fitting procedure (and 
so from the error bars of the data) and from the variation 
of the second-multiplet mass My 2 inside the considered 
region (6). Our fit to the NA60 data has a y 2 per degree 
of freedom of 1.5. 

Apparently we obtain a rather satisfying fit to the NA60 
data except for the last two or three data points at the 
largest values of the dimuon mass. At present, none of 
the hadron-theory approaches to this w transition form 
factor [31, 47, 39, 48, 41] is able to understand these last 
data points. Recently it has been suggested in [49] us¬ 
ing a dispersive calculation and high-energy constraints 
that these data points might be incompatible with QCD. 
Clearly it would be highly desirable to obtain additional 
data for this ui transition form factor, in particular from 
experiments where the complete final state n°£ + £~ can be 
reconstructed and not only the dilepton; see [28] for more 
details. 



0.1 0.2 0.3 0.4 0.5 0.6 

V?[GeV] 

Figure 2: The normalized u-n transition form factor as a func¬ 
tion of the (time-like) virtuality. Various theoretical calculations 
are displayed together with NA60 data obtained from the decay 
uj —► [28]. See the main text for more details. 


In the following we will stick to our result for k as given 
in (31) and obtained from a fit to the full range of NA60 
data since it is the best that one can do in the present 
situation. However, we briefly discuss two alternatives: If 
we performed an alternative fit to the NA60 data rejecting 
the last three data points, we would get k = 19 ± 2 , which 
is fairly compatible with our full fit. The y 2 per degree of 
freedom would reduce to 0.8. If we rejected the NA60 data 
altogether and regarded the dispersive calculation of [39] 
as the “truth”, we would find n = 13.1 ± 0.5. The result of 
this last fit is also shown in Fig. 2 as the dark gray band. 
If we return to the first panel of Fig. 1 we observe that 
one obtains the same qualitative features with such values 
of k. Nevertheless, we will not use these results further 
on and will take into account only the result (31) of the 
all-data fit. 

As for the pion transition form factor we shall compare 
our result to similar approaches from the literature. The 
form factor from [31], 

- TI „ M 2 ± <7 2 

^ y (g 2 ) = M V 2_ q2 , (32) 

labeled by “TL” in Fig. 2, lies within our uncertainty band 
in spite of the fact that the derivation is based on a very 
different approach. The results for VMD and LMD can be 
obtained from (20) and (21), respectively, using (28) and 



































Figure 3: Feynman diagrams for the cj —>■ 7r°7 and to —> 7 
decays at leading order in QED. The shaded blob corresponds to the 
■^ 7 r°cj (*) (^ 2 ) f° rm factor and £ denotes e or ^ l. 


(30). For instance, the VMD expression is 


fVMD 
J 7T°U)V 


( 9 2 ) 


My r - q 2 ■ 


(33) 


As can be seen in Fig. 2 the results for both VMD and 
LMD deviate significantly from the data and therefore 
from our approach. They also deviate from the results 
of the dispersive calculation of [39], which is shown as the 
black band in Fig. 2. As already noted we have performed 
an alternative fit of our THS expression to the dispersive 
result — the dark gray band in Fig. 2. It should be noted 
that our rational function cannot fit the cusp structure 
in the low-energy region yfq* £ [0.3, 0.4] that emerges in 
the dispersive calculation from a cross-channel inelasticity; 
see [39] for a detailed discussion. We repeat our statement 
that additional data on this transition form factor would 
be highly welcome. 

To further explore the validity of the THS scheme we 
study selected decay channels of the ui meson. It appears 
to be convenient to introduce the matrix element 


MZoM = 


ZuiZir 


e^ ap Paqp 


lim (p 2 - Afy )r 2 n(r 2 ;p 2 , q 2 ) . 

r z —> 0 


(34) 


Using previous definitions (10) and (28) this can be re¬ 
duced to 


(p> q) = e^°Mq*V vaf) P a q 0 ■ (35) 

In (34) the factors 2/3 and 1/3 emerge from the respec¬ 
tive overlap of the u meson with the singlet and the 8th 
component of the octet current, i.e. the factors come from 
(3). These factors actually sum to unity, which is not sur¬ 
prising from a flavor SU(2) point of view where the a; is a 
pure singlet. 

We start with the prediction of the decay width of the 
ui —>• 7r°7 process, which is depicted in Fig. 3a. The matrix 
element 


o 7 = M^ 0 {p, q)e ll (p ., A w )e*(<f, A 7 ) 


(36) 


does not depend on re as already mentioned. This provides 
us with a pure prediction. After simple manipulations we 
find 


r(w 


0 olMI 
V 7) = —rn K 7 t° 


24 

o?Z 2 LF, 


wv(0)| ^ 

°uiV ( 0 ) | 2 


1 - 


162 r(, 


UJ 




Ml 
M*. 
M 2x3 


1 - 


Ml 


(37) 


where we have used (26) to obtain the last expression. 
Taking N c = 3, F = (92.22 ± 0.06) MeV [50], = 

134.98 MeV, M Vl = M p = (775.26 ± 0.25) MeV [22] and 
My 2 in the range (6) we obtain the values listed in Tab. 1. 
We see that THS is within its uncertainties compatible 



I VMD 

LMD 

THS 

experiment 

r(u - 

4 7T°7) [MeV] | 0.68 

0.45 

0.63 ± 0.04 

0.70 ± 0.03 


Table 1: Decay width of the process uj —> 7r°7 as obtained from 
different approaches and from data [22]. We estimate an uncertainty 
only for our THS case. 


with the experimental data. In this case VMD proves to 
be phenomenologically successful, too. The value of the 
LMD approach seems to be off even taking into account 
the 40% rule-of-thumb uncertainty used in [17]. 

Moving on to the oj —> tt°£ + £~ decay, which is depicted 
in Fig. 3b, the matrix element can be written in the form 

(_*) _ 

= M l l' n 0 (p,q)e IJ .(p,\ u ,)Z- r u(q 1 )(-ie'Y„)v(q2). 

(38) 

The decay width can then be expressed as a form-factor- 
dependent integral over the dilepton invariant mass, 




r(w ->• n°£ + r) = 


TW? 2 )] 2 


72t tMI 


4 m? 


(39) 


X 4/1 - 


4m? 


1 + 


2 m 2 


Ai(M 2 ,M 2 , 9 2 )dg 2 , 


where A denotes the Kallen triangle function defined as 
A (a, 6, c) = a 2 + b 2 + c 2 — 2 ab — 2 ac — 2 be. (40) 


For THS the result for this decay width depends on the 
parameter re. We use the fitted value (31) and for the other 
input quantities the ranges specified after (37). 

It is common practice to normalize decay widths that 
involve dileptons to the corresponding decay widths in¬ 
volving photons [30]. This leads to the branching ratios 
listed in Tab. 2. Obviously, for the electron case all the 
results of the considered approaches lie within the experi¬ 
mental uncertainty. For the muon case, however, only the 
THS result explains the experimental value. This should 
not come as a surprise given that in Fig. 2 the LMD and 
VMD curves lie much lower than the THS band. 


q 2 —0 
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VMD 

LMD 

THS 

experiment 

B(^->,O e +e-) .. i A 3 
x 10 

9.1 

8.9 

9.6 ± 0.1 

9.3 ± 1.0 

) x 10 3 

0.91 

0.82 

1.33 ± 0.10 

1.6 ± 0.5 


Table 2: Branching ratios of the decays u 7r°e"*"e — and cj —> 
7rnormalized to the branching ratio of the decay lj —> ir 0r y as 
obtained from different approaches and from data [22]. We estimate 
an uncertainty only for our THS case. 


4. The process 7r° —> e + e 



Figure 4: Leading-order contribution to the 7r° —> e+e — process in 
the QED expansion and its representation in terms of the leading 
order in chiral perturbation theory. The shaded blob corresponds to 
the doubly virtual pion transition form factor (l 2 , (r — l) 2 )- 


In the following two sections we turn our attention to 
the decays of the neutral pion. Namely we will discuss the 
rare decay 7 r° —> e + e~, to which this section is devoted, 
and later on also the Dalitz decay 7r° —> e + e~ r y. 

The decay 7r° —> e + e _ and the radiative corrections 
connected to this process have been extensively studied 
in [6, 7]. On the left-hand side of the graphical equation 
in Fig. 4 the leading order of the QED expansion of the 
considered process is depicted. The doubly virtual pion 
transition form factor J r 7r o 7 . 7 . (/ 2 , (r — l) 2 ) is represented 
in the figure as the shaded blob. Here r denotes the pion 
momentum and l a loop momentum. This transition form 
factor plays here an essential role, since it serves as an 
effective ultraviolet (UV) cut-off due to its l/l 2 asymp¬ 
totics governed by the OPE. The loop integral over d 4 l is 
therefore convergent. 

On account of Lorentz and parity symmetry the on-shell 
matrix clement of the 7r° — > e + e - process can be written in 
terms of just one pseudoscalar form factor in the following 
form: 


iM^o^e+e- = P 7T o^ e+e - [u(q 1 )'y 5 v(q 2 )} . (41) 


Subsequently, the decay width reads 


r(7r° 


e + e 


) = 


2M 2 \P 7r o 


— >-e+e - I 


107 tMtt 


II- 


4m 2 

M2 


(42) 


Taking into account only the leading order (LO) in the 
QED expansion, i.e. the left-hand side of the graphical 
equation in Fig. 4, we find (within the dimensional reduc¬ 
tion scheme [51, 52]) 


ie 4 m e r d 4 l P^o^ T (p 2 ,q 2 )X{M 2 ,p 2 ,q 2 ) 
M 2 J (27t) 4 p 2 q 2 (l 2 — ml) 

(43) 


Here p = l — q\ and q = l + q 2 , where qi and q 2 are the 
lepton momenta. Using dimensional regularization and 
Passarino-Veltman reduction [53] the explicit result of the 
loop integration in terms of scalar one-loop integrals is 
given in Appendix A for various form factors. 

On the right-hand side of Fig. 4 we can see how the pre¬ 
viously discussed expression is represented in leading order 
of chiral perturbation theory [16, 17]. Here, the constant 


Tr ° 7 * 7 *( 0 , 0 ) 


N c 

12t r 2 F 


(44) 


emerging from the chiral anomaly is used instead of the 
full transition form factor; see (16). This leads to a diver¬ 
gent integral [15] and it is thus clear that a counter-term 
is needed. If the loop is renormalized at scale p, the fi¬ 
nite part of such a counter-term is denoted by x^(aO- It 
corresponds to the high-energy behavior of the complete 
transition form factor. Therefore, concerning the pion de¬ 
cay, we can characterize each transition form factor by the 
corresponding value of X^Hp)- 

In fact, using (44) instead of the full form factor in (43) 
together with the counter-term Lagrangian from [16, 17] 
one obtains 


e£L+ e - = -2a' J m e 7vo r7 .{0,0) { Q + ^ log (A f 




Li 2 (z) — Li 2 ( — 


- J7rlog(— z) 


(45) 


In the above formula Li 2 is the dilogarithm, j3 e = 
-\/l — Aml/M 2 , z = — (1— /3 e )/(l+/3 e ) and p represents the 
scale at which the loop integral is effectively cut off; cf. [17]. 
This effective approach with the constant form factor lead¬ 
ing to formula (45) can be conveniently used to match 
with the calculation (43) and any full, i.e. momentum- 
dependent form factor. Thus, for various approaches, the 
value of the associated effective parameter x^ (p) can be 
extracted. In other words, the left-hand side of (45) can 
be substituted by the expressions for P^P_+ e+e - stated in 
Appendix A and X^(p) I s subsequently determined. Fol¬ 
lowing common practice [11, 46, 16, 17] we have chosen 
p = 770 MeV ~ M p to set the renormalization scale. 

An approximative way to get the low-energy constant 
x (r \p) up to corrections proportional to m e /M n and 
M n /p has been presented in [4]. In our notation this reads 


x (r V) 


5 

4 


d t log (t) 


d J r 7r o 7 * 7 . (— tp 2 , — tp 2 ) 

dt J> 7 , 7 ,(0,0) 


(46) 

Using this approach we get very simple formulae for the 
VMD model 


Xvmd(m) - Y “ \ lo § 
as well as for the LMD case (cf. [17]) 

(r) / \ 11 3 M^ 

Xlmd(m) - -j - 2 l0g ~iP~ ’ 


K 

i, 2 


127 t 2 F 2 

iVcM^‘ 


(47) 


10 


(48) 























For the THS form factor, where the additional mass scale 
My 2 and the parameter n come into play, the calculation 
according to (46) yields 


(l) f \ 5 3 M V 1 , 3 My My 

Xths(^) - i ~ 2 1 g + o 717!>-775 


^k 

Mk 

^k 

Mk 


Ml 14 

—-(27 tF) 2 

Ml_ N r y ! 


l v 2 
2 M* ± 


Ivl v 2 

1 


k 2 


7 


(4? tF) 2 


« K, M k 

iV c (4^F) 4 


My 2 ~ My ± N c (M* a -Mk )J 

M 2 2 + ^ 


, A4 a 
l0§ ^ 


2 (M^-M^ 


■log 




- 1 


(49) 


the photon-energy cut I £ 7 < M (1 — x cut )/2 in the rest 
frame of the decaying pion, with x cut = 0.95. Soon after, 
the disagreement of (50) with a theoretical calculation was 
found [4], 

To predict the quantity (50) we proceed as follows. Tak¬ 
ing into account NLO QED corrections [ 6 , 7] we can cal¬ 
culate the branching ratio of the inclusive process 7r° —>• 
e+e" ( 7 ), i.e. of a process where we allow bremsstrahlung 
photons to appear in the final state. Denoting these QED 
corrections by S this can be written as 


B(tt° -i e + e ( 7 ), x > ar cut ) 

_ r LO (7T° e+e~)[l + <5(s cut )] 
r L o( 7r o —> 77) 


x B(tt° —1 77 ). 


(51) 


The numerical values for the exact as well as for the 
approximate results are listed in Tab. 3. In the first row the 


In addition, the width of the pion decay into two photons 
has the form 


(p = 770 MeV) 

VMD 

LMD 

THS 

exact result 

2.87 

2.29 

2.2 ± 0.2 

approx, from (46) 

2.72 

2.2 ± 0.9 

2.0 ± 0.2 


Table 3: Values of the effective parameter (770 MeV) for vari¬ 
ous pion transition form factors. In general we provide uncertainty 
estimates only for THS. The uncertainty for LMD stems from [17]. 

values gained from the matching of (43) to (45) are listed. 
The second row contains the results of the approximative 
formulae (47), (48) and (49). In the LMD case, the value 
according to (48) would be 2.16, but here the listed result 
has been taken from [17] where also a 40 % uncertainty has 
been estimated. Needless to say, the used numerical inputs 
were the same as in the whole text except for one thing. 
In the case of the pion decays, we use for the ground-state 
vector-meson multiplet mass My 1 the average of the p and 
oj meson masses, i.e. My x = (779 ± 4) MeV. 

We can see that the approximative formula (46) is in¬ 
deed reasonable. The numerical values are close to the case 
of the exact calculation, where electron and pion masses 
are not neglected. To compare briefly the values for the 
various approaches, we see in Tab. 3 that the THS mean 
value is close to LMD but incompatible with the VMD 
value. Of course, a meaningful comparison would require 
an uncertainty estimate for the VMD case, which is not 
available. Considering the large uncertainty of LMD ac¬ 
cording to [17], THS and even VMD lie safely within the 
LMD band. In turn the central value of LMD lies within 
the THS band. 

Finally, we can utilize the above formulae and pro¬ 
vide branching ratios to be compared to experiment. The 
presently most precise measurement was performed by the 
KTeV collaboration at Fermilab with the result [3] 

B KTeY (Tr° e+e" ( 7 ), x > 0.95) 

* (50) 

= (6.44 ± 0.25 ± 0.22) x 10 " 8 , 

where x is the normalized lepton-pair invariant mass, i.e. 
x = (?i + < 72 ) 2 /M%. This condition can be translated to 


r LO ( 7 r° —»• 77 ) 


1 1 e 4 M 4 

2 16?rM T 2 


|J r 7r o 7 . 7 .(0,0)| 2 


(52) 


and B(tt° ->■ 77 ) = (98.823 ± 0.034) % [22], The calcula¬ 
tion of the two-loop virtual radiative corrections together 
with the bremsstrahlung correction in the soft-photon ap¬ 
proximation gives the result <5(0.95) = (—5.8 ± 0.2) %. It 
has been shown in [7] that the soft-photon approximation 
is a valid approach for the value x cu t=0.95 used by the 
KTeV experiment. Note that this value of the correction 
S was obtained in [ 6 ] in a model-independent way using 
y( r ) (770 MeV) = 2.2 ± 0.9. However, this result for S de¬ 
pends negligibly on the range of considered values of 
shown in Tab. 3. Hence we consider the result for <5 to be 
valid for our case. 

With this NLO QED input we find the values shown in 
Tab. 4. Here we see that the central value in the THS 


| VMD LMD 

THS 

KTeV 

B(tt° —»• e+e“( 7 ) , x > 0.95) | 5.96 5.8(3) 

5.76(7) 

6.44(33) 


Table 4: Branching ratio of the inclusive process 7 r° —> e+e — ( 7 ) at 
NLO of QED for various models for the pion transition form factor. 
The listed values are to be multiplied by a factor 10 -8 . The KTeV 
value is based on the result stated in [3]. 

case is, of course, compatible with the LMD central value, 
since the same holds for the values x^(/-0; see Tab. 3. 
The uncertainty of the LMD approach according to [17] 
makes it compatible with the KTeV result on the level of 
1 cr. Considering VMD, where we do not have an estimate 
of the theoretical uncertainty, we are then 1.5 a away from 
the experimental value. Finally, the THS approach is 2 <7 
off. As compared to LMD this is related to the fact that 
the estimate for our theoretical uncertainty is on the level 
of only 1 % concerning this particular branching ratio. We 
recall that this estimate is merely based on two sources: 
first, the fit for k, which might be improved in the fu¬ 
ture by data on the doubly virtual pion transition form 
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factor; second, a large variation in My 2 . By this varia¬ 
tion we estimate the uncertainty imposed by a truncation 
of the infinite tower of resonances. Strictly speaking our 
model intrinsic uncertainty estimate does not account for 
the two initial assumptions: the chiral and the large-A/, 
limit. We expect that to some extent we have accounted 
for the uncertainties of the large-_/V c approximation by fo¬ 
cusing on kinematical situations where the (large--W sup¬ 
pressed) widths of the resonances do not matter so much 
and by varying the mass of the second vector-meson mul- 
tiplet in a relatively large range. The uncertainties caused 
by the chiral limit should be on the order of M 2 / A 2 where 
A x denotes the scale of chiral symmetry breaking. Within 
chiral perturbation theory one uses the scale where new de¬ 
grees of freedom come into play or where the loops become 
as important as the tree-level contributions. Roughly this 
implies A x < M p ,47rF, i.e. an uncertainty of about 3%. 
Adding errors in squares this is clearly negligible as com¬ 
pared to the 10 % uncertainty documented in Tab. 3 for 
THS. Based on these considerations we regard our uncer¬ 
tainty estimate in Tab. 4 as reasonable — assuming that 
our determination of k from the uj-tt transition is correct. 

Consequently, if the (central value of the) KTeV re¬ 
sult will be confirmed by future experiments and mea¬ 
sured with higher precision such that the discrepancy 
would reach several er’s, then the following two scenar¬ 
ios are conceivable: Either there are some aspects of the 
THS approach which are not well-suited for the rare de¬ 
cay 7r° —> e + e ~, or it should be seriously considered that 
physics beyond the Standard Model influences the rare 
pion decay to a significant extent. However, under the 
present circumstances we consider the current discrepancy 
to be inconclusive. We will come back to this discrepancy 
in the next section. 


5. Singly virtual pion transition form factor 
and the Dalitz decay i r° —> e + e~j 


With particular models at hand, we can also explore the 
properties of the singly virtual pion transition form factor 
Jq r o 7 * 7 . (0 ,q 2 )- A slope a T and a curvature b w of the form 
factor are defined in terms of the Taylor expansion in the 
invariant mass of the vector current [40] 


j~ 7T°7*7* (0; Q ) 

J' 7r o 7 . 7 »(0,0) 


~ 1 


7 Ml 



2 


(53) 


By means of the first and second derivative it is straight¬ 
forward to find the following expressions for the simplest 
form factors 


VMD 

a 7T 


Ml 
Ml ’ 


M 2 

7 VMD _ lvl ir VMD 

w “ Ml a * ’ 


(54) 


a 


LMD 

7T 


Ml r 4 tt 2 F 2 ' 

MfX ~ NlMfJ ’ 


iLMD 

^7T 


Ml lmd 

Ml * 


(55) 


as well as for the THS form factor with its two vector- 
multiplet mass scales 


THS 


6 THS = 


Ml 
Mh 


1 + 


Mi 


V i 


Ml 


Ml 

Ml 


M, 


2 1 


Ml 


24t t 2 F 2 


THS 


Mi 


(56) 


Ml Ml 


Numerical results are shown in Tab. 5. We see that THS 



VMD 

LMD 

THS 

dispers. [40] 

exp. [22] 

slope a n 
curvature b^ 

30.0 

0.90 

24.5 

0.74 

29.2(4) 

0.87(2) 

30.7(6) 

1.10(2) 

32(4) 


Table 5: Slope and curvature of the singly virtual pion transition 
form factor evaluated for various approaches. The listed values are 
to be multiplied by a factor of 10 —3 . 

is compatible with the experimental value and that it is 
numerically consistent with VMD. Similar to some pre¬ 
vious cases, the value predicted by LMD is unsatisfying. 
Nevertheless, none of the models (VMD, LMD, THS) fully 
agrees with the results of the dispersive calculation [40]. 
The latter is capable of producing very reliable results at 
low energies. We note in passing that the apparent agree¬ 
ment of VMD with the dispersive result for the slope is 
somewhat accidental. Actually the discrepancy grows, if 
the <f> meson is included in the analysis [40]. 

The disagreement between THS and the dispersive re¬ 
sults brings us back to the discussion of the uncertainties 
inherent to THS and to our attempt to quantify these un¬ 
certainties in a reasonable way. Of course, for a given 
model one can only estimate the model intrinsic uncer¬ 
tainties, which in turn inherit a model dependence. Our 
numerical values for the uncertainty estimates of the THS 
results emerge dominantly from the uncertainty in the de¬ 
termination of the parameter n and from a variation in 
My 2 . In contrast to most of the results previously pre¬ 
sented, the singly virtual pion transition form factor does 
not depend on the parameter k. This can be most eas¬ 
ily seen from (56), but also from (15), (17). In addition, 
the dependence on My 2 is minor for the low-energy quan¬ 
tities determined in (56). In fact, in (56) the mass My 2 
essentially shows up in the numerically small combination 
AM 2 /Ml with AM 2 = Ml - 24t : 2 F 2 /N C . The small¬ 
ness of AM 2 might be seen as an incarnation of the KSFR 
relations; see, e.g., [54] and references therein. In turn, 
this implies that our model intrinsic uncertainty determi¬ 
nation might underestimate the real uncertainty if the de¬ 
pendence of our results on My 2 and/or k is accidentally 
small. While this is the case for the singly virtual pion 
transition form factor, it is certainly not true for the dou¬ 
bly virtual one; see also Fig. 1 with the broad uncertainty 
band in the first panel and the corresponding nearly in¬ 
visible spread in the second panel. Obviously, one has to 
take our uncertainty estimates with a grain of salt. 
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(a) (b) (c) 

Figure 5: Leading order diagram of the Dalitz decay n° —► e+e — 7 
in the QED expansion (a) and the one-photon irreducible contribu¬ 
tions to the NLO virtual radiative correction (b) and (c). Note that 
“cross” in (b) accounts for a diagram with a photon emitted from 
the outgoing positron line. Diagrams (b) and (c) also serve as a 
bremsstrahlung contribution to the 7r° —> e+e — decay. 


As a closure of the phenomenological part of this work, 
we inspect the theoretical predictions for the decay width 
of the Dalitz decay 7 r° —» e + e~ r y. The leading-order QED 
contribution is depicted in Fig. 5a. Our aim is also to ad¬ 
dress the NLO radiative corrections. For this purpose we 
use the approach documented in [42], which has recently 
reviewed and extended the classical work of [55]. Hence, 
together with the bremsstrahlung beyond the soft-photon 
approximation we use in the following calculations the vir¬ 
tual radiative corrections including also the one-photon ir¬ 
reducible (I 7 IR) contribution. The diagrams connected 
with the latter part are schematically shown in Fig. 5b 
and 5c. The mentioned correction then emerges when 
one considers the interference term with the LO diagram 
in Fig. 5a. 

The decay width of the Dalitz decay at NLO can be 
written in the following form 

r( 7 r° - 4 - e+e“ 7 ( 7 ), x > x cut ) 

TO n , (57) 

= [l + S D (x> x cut )] r u (tt u ->■ e + e 7 , x > x cut ), 


where we have denoted the sum of all the NLO radiative 
corrections to the integrated decay width as Sd and where 


r LO ( 7 T° -»• e + e 7 , x > x cut ) 

= ^ — ^ r LO (7T° —> 77 ) f dx 

“ & cut 


X 


8 Pe{x) (1 - x ) 3 / 2 ml \ 

3 4x \ xM%) 


•7vy 7 «(0 ,xMp 

'7y°7*7* (O 5 0) 


2 


(58) 


In the above formula we have introduced /3 e (x) = 

i/l — £rri*/(xM%). Note that contrary to the last section 
x cu t cannot be interpreted as a cut on a single photon 
energy, since the additional bremsstrahlung photon ap¬ 
pears. There are also no restrictions on the energy of the 
bremsstrahlung photon in this decay width. 

The numerical results for the considered approaches are 
compared in Tab. 6 . In the first row the NLO correc¬ 
tion is presented. Naively one might get from the defi¬ 
nitions (57) and (58) the impression that the pion tran¬ 
sition form factor does not enter Sd- This is, however, 



WZW 

VMD 

LMD 

THS 

NLO correction <5j) 

8.338 

8.299 

8.307 

8.3008(3) 

r(7T°->-e+e - 'y(-7)) 

r LU (7rO-).-Y-Y) 

11.9503 

11.9735 

11.9692 

11.9729(4) 

r(7T 0 ->-e + e - 'y('y) , x>0 .23191 

r(7r u — >e + e — 7 ( 7 )) 

31.03 

31.63 

31.52 

31.61(1) 


Table 6 : Theoretical prediction for the Dalitz decay branching ratios 
at NLO of QED with and without a cut as evaluated for various 
approaches. “WZW” describes the case where the pion transition 
form factor is replaced by its low-energy limit. The listed values are 
to be multiplied by a factor of 10 —3 . 


not the case. The NLO correction involves, of course, the 
four-body phase-space integrations and the pion transition 
form factor cannot be scaled out. This is due to the in¬ 
tegration over the x variable, which is tacitly performed 
also on the left-hand side of (57). On the other hand, 
taking into account the correction to the differential de¬ 
cay width, we would get an expression independent of the 
form factor. This statement then holds only when leaving 
out the I 7 IR part, where the form-factor-dependent loop 
integration is performed. Nonetheless, concerning model 
uncertainties, it is encouraging that numerically the NLO 
correction is to a very large extent independent of the used 
pion transition form factor. Using the respective value of 
the NLO correction, the second row in Tab. 6 shows the 
corresponding branching ratio if no cut on additional pho¬ 
tons is applied, i.e. the branching ratio of the Dalitz decay 
for the whole kinematic region x € [Arri^/M*, 1]. The re¬ 
sults of this row should be compared to the experimental 
value = (11.88±0.35) x 10 -3 , which has been 

calculated using inputs from [22]. The agreement with the 
values of the second row in Tab. 6 is very satisfying. 

In the previous section we have stated the value (50) as 
the outcome of the KTeV analysis. However, what KTeV 
has determined at first place is 3 


r( 7 T° —> e + e ( 7 ), x > 0.95) 
r( 7 T° —> e+e _ 7 ( 7 ), x > 0.2319) 

= (1.685 i 0.064 ± 0.027) x 1CT 4 . 


(59) 


We note in passing that for the quantity displayed in the 
last row of Tab. 6 the corresponding value used by KTeV 
to obtain finally the result (50) is 3.19% [3]. This value is 
close but not equal to any of the results in the last row of 
Tab. 6 . 

Having the radiative corrections determined from the¬ 
ory suggests to calculate directly the ratio (59) instead of 
the derived quantity (50) and to compare the calculation 
to the experimental result. However, the result (59) is not 
a direct measurement either, but is, of course, based on 
an interplay of measurements, Monte Carlo simulations, 
acceptance corrections, etc. [3, 56]. By the time the KTeV 


3 Note that we use here the more precise value x cu t = 0.2319 

from [56] in contrast to 0.232 as stated in [3]. This value for Xcut 
corresponds to (65MeV/M*-) 2 . 
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measurements were analyzed not all QED radiative correc¬ 
tions were available. To evaluate the decay width of the 
Dalitz decay the radiative corrections from [55] have been 
used. In particular the I 7 IR contributions have not been 
included in [55] . In addition, in the Monte Carlo generator 
of KTeV the radiative corrections for the e + e _ decay have 
been adopted from [57]; see [56]. However, in [57] some 
approximations have been used that have been proven to 
be misleading after the exact calculation has been per¬ 
formed [ 6 ] . Therefore it would be somewhat misleading to 
compare a theory result with all radiative corrections in¬ 
cluded to an experimental result where only part of these 
corrections have been taken into account. On the other 
hand, it might be illuminating to provide a rough esti¬ 
mate of the impact of these differences. Before presenting 
such an estimate we would like to repeat the take-home 
message from the first section: All NLO QED radiative 
corrections are now available [ 6 , 7, 42] and can be taken 
into account in future analyses of data on the pion decays 
7 r° —>• e + e _ and 7r° —» e + e _ 7 - 

In spite of the fact that not all radiative corrections 
have been taken into account in the analysis that led to 
the result (59) we do not feel legitimated to modify this 
experimental result. Instead we will calculate the ratio of 
(59) within THS 4 , but assign an uncertainty to it, which is 
related to the neglect of the I 7 IR contributions . 5 It should 
be stressed, however, that this is not an entirely valid ap¬ 
proach either. After all, as soon as one can calculate the 
NLO QED radiative corrections, they do not constitute 
an uncertainty of a result but rather a well-defined shift of 
the result. Nonetheless, we proceed in the described way 
to obtain a rough estimate of the uncertainty caused by 
neglecting part of the radiative corrections. 

The quantity of our interest can be generally expressed 
and expanded as 


A±tr(A) 

B±a{B) 


A 

B 


1 + T a ( B ) 

A + B 


(60) 


In our case we have A = r(7r° — > e + e~( 7), x > 0.95) and 
B = r(7T° —> e + e _ 7 ( 7 ), x > 0.2319), which we calculate 
up to NLO using the complete set of corrections. The un¬ 
certainties are then given by the I 7 IR contributions to the 
NLO virtual radiative corrections, i.e. by the interference 
of the diagrams shown in Fig. 5. Hence 

cr(A) _ |r^R(7r° -7- e+e“ 7 , x > 0.95)| 

A r(7T° -7- e+e _ ( 7 ), x > 0.95) 

|£d 7IR (x > 0.95)| r LO (7r° -7 e + e“ 7 , x > 0.95) 

“ [1 +<5(0.95)] r LO (7T° -7 e+e-) ’ 

(61) 


4 According to Tab. 6 the results from the other approaches are 
comparable. 

5 The modification induced by using [57] in the Monte Carlo gen¬ 
erator is even harder to assess. Therefore we concentrate solely on 
the I 7 IR contributions. 


which numerically yields &(A)/A ~ 0.4 %. For the second 
term we get 

a{B) |r™j°( 7 r° -7 e+e- 7 , x > 0.2319)] 

"~B"~ " r(7T° -7 e+e- 7 ( 7 ), £>0.2319) 

|(5^' yIR '(cc > 0.2319)] 

“ [1 + Sd(x > 0.2319)] ’ 

for which we find cr(B)/B ~ 0.5%. We can thus con¬ 
clude that the total uncertainty would be at the level of 
1 %. For completeness we list the values that we have 
used to get the above stated numbers. The input for (61) 
is (5p 7 lR (a; > 0.95) = —8.93(1) % for the I 7 IR contribution 
to the NLO virtual radiative correction for the decay width 
of the Dalitz decay. For the ratio of the LO decay widths 
of the Dalitz decay with x > 0.95 and the 7r° -7 e + e“ 
process we have -— ( p L o ^°%2+e ~)°' 95) = 4.31(5)%. In 
(62) we have utilized for the overall NLO correction to the 
Dalitz decay Sd(x > 0.2319) = —3.189(2)% and finally 
<5 d 7IR (:e > 0.2319) = -0.5038(3) %. In general, we can see 
that the radiative corrections are getting more important 
with higher values of x cu t- Taking values from Tab. 4 and 
Tab. 6 and the previously estimated relative uncertainty 
of 1 % (we state this separately) we find 


r( 7 r° —» e + e ( 7 ), x > 0.95) 
r( 7 r° —> e+e _ 7 ( 7 ), x > 0.2319) 


theo+uncert. 


= (1.54 ±0.02 ±0.02) x 10" 4 . 


(63) 


This result is to be compared with the KTeV value (59), 
yielding a discrepancy of approximately 1.5 a. We stress 
one more time that the result (63) should be read with 
caution. 


6. Outlook 

Several directions are conceivable how the THS scheme 
might be extended or combined with other approaches. 
Most straightforward would be to extend THS to other 
correlators; see also [5]. Eventually one might get cross 
correlations and in that way more information, e.g., by 
fitting to additional sets of data. 

In the present work we have merely used the PVV cor¬ 
relator as an intermediate step to obtain the nVV corre¬ 
lator. Of course, the PVV correlator is also interesting in 
its own right. For instance, as it has been worked out in 
[5] , the low-energy expansion of this correlator determines 
some of the low-energy constants of the chiral Lagrangian 
at order p e in the anomalous sector [58], i.e. the correc¬ 
tions to the Wess-Zumino-Witten Lagrangian. Having 
determined some of the parameters of the PVV correla¬ 
tor within THS by demanding high-energy constraints one 
could provide estimates and/or cross-correlations for these 
low-energy constants following the procedure outlined in 
[5], 
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As already spelled out in Section 1 we refrained from 
using high-energy constraints for quantities like the tujjV 
correlator since the vector meson might probe already too 
many details of the intermediate-energy region (above 1 
GeV) where we have approximated the tower of infinitely 
many states by one effective state per channel. Extending 
the scheme from THS to three hadrons per channel, albeit 
introducing more parameters in the first place, might allow 
to use quark-scaling relations along the lines of [18, 19, 20] 
for 7 tujV and similar correlators which in turn could help 
to keep the number of free parameters manageable. In 
the same spirit one can use subleading orders in the high- 
energy expansion in terms of QCD condensates [43]. 

In the present work we have used the chiral limit to 
construct the PVV correlator. In line with this approxi¬ 
mation we have focused on the nVV and ttloV correlator 
where no strange quarks are involved. For a reasonable 
extension from the pion to the whole pseudoscalar multi- 
plet (a nonet in the large- N c limit) one has to go beyond 
the chiral limit. The prospect of such an endeavor would 
be the possibility to tackle the transition form factors of 
the rj and rf mesons and their rare decays into electron 
and positron or muon and anti-muon. The simplest way 
to go beyond the chiral limit would be to utilize a per¬ 
tinent Lagrangian. In extensions of chiral perturbation 
theory [8, 11, 12] one frequently uses a Lagrangian where 
the vector mesons are represented by antisymmetric tensor 
fields [11, 46, 37, 31, 32, 33, 35, 36]. As discussed in [5] a 
given resonance Lagrangian of formal chiral order p e might 
not be capable of satisfying all desired OPE constraints. 
Indeed, as shown in [37] this problem appears at least if 
only one vector-meson multiplet is used for the PVV corre¬ 
lator. Concerning the vector instead of the antisymmetric 
tensor representation the same negative result has been 
reported in [5]. It might appear that the first-order ap¬ 
proach suggested in [59] offers more flexibility here. Alter¬ 
natively or in addition one might consider the inclusion of 
resonance terms of chiral order p s . The following consid¬ 
eration might demonstrate that one can expect to create 
additional terms that would not spoil the OPE right away: 
If one normalizes the PVV correlator to the correspond¬ 
ing chiral-anomaly expression, a p 6 resonance Lagrangian 
would produce “corrections” that schematically are given 
by 


£p 6 ,R —•► 



with n < 3 . 


(64) 


We have replaced all resonance masses by Mr and all mo¬ 
menta by Q. This is sufficient to discuss the high- and 
low-energy limits. The factor of Q 2 in the numerator of 
(64) emerges from the fact that the resonance Lagrangian 
is assumed to be of chiral order p 6 while the chiral anomaly 
is of order p 4 . The limitation n < 3 comes from the fact 
that a three-point correlator is considered. At high ener¬ 
gies a single term with n = 1 definitely spoils the OPE. 
Coupling constants between various contributions must be 


adjusted such that cancellations appear. Of course, this is 
the reason why the OPE is capable to tell something about 
the resonance parameters. A resonance Lagrangian of chi¬ 
ral order p s produces terms with the schematic structure 


C 


P S ,R 



with n < 3 . 


(65) 


Here the terms with n = 1,2 spoil the OPE and must 
cancel. The terms with n — 3, however, produce struc¬ 
tures that could be compatible with the OPE but are 
structurally different from the ones in (64). Thus the in¬ 
clusion of resonance terms of chiral order p s might help 
to obtain an agreement between a Lagrangian approach 
and the OPE. Needless to say that the systematic con¬ 
struction of a p s Lagrangian would be extremely tedious; 
see [37] and references therein. In practice, according to 
[21], leading-order OPE constraints 6 together with a singly 
virtual pion transition form factor that vanishes at large 
momenta can be achieved by a resonance Lagrangian with 
two vector-meson multiplets in the antisymmetric tensor 
representation. Note, however, that the concrete B-L limit 
(14) has not been used in [21]. Nonetheless, the results of 
[21] suggest that THS can be reproduced from a resonance 
Lagrangian of order p 6 in the chiral counting. 

While leaving the chiral limit might be systematized by 
a Lagrangian approach, our second basic assumption, the 
large-iV c limit is much harder to remedy. On purpose we 
stuck to narrow states (7r, uj) and to the space-like or low- 
energy time-like region where the effects of inelasticities 
are minor. Nonetheless this induces uncertainties, which 
are hard to assess on a quantitative level. At low ener¬ 
gies an excellent way to deal with inelasticities is to use 
them in ones favor by a dispersive setup [39, 40, 49]. A 
powerful framework emerges when combined with chiral 
low-energy constraints, QCD high-energy constraints, and 
precise data on the scattering amplitudes that come into 
play via the inelasticities. In practice the scheme becomes 
intractable at intermediate energies (roughly above 1 GeV) 
where more and more many-particle channels open up. 
Even if one studies low-energy quantities one might need 
some information from the intermediate-energy regime if 
the QCD high-energy limit is reached very late. This is 
exactly what we saw for the THS result of the 7rW cor¬ 
relator. This suggests to study the impact of the THS 
results on a dispersive calculation by using THS in the 
intermediate- and high-energy space-likc regime where its 
results should be most trustworthy. 

In the present work we have determined the 7 rVV corre¬ 
lator with a focus on the rare pion decay into electron and 
positron. In turn this decay might offer some window to 
observe low-energy traces of physics beyond the Standard 
Model. However, the nVV correlator is also interesting 
because it contains information about the intrinsic struc¬ 
ture of the pion. Hence it could be illuminating to figure 


®Note that “chiral order” refers to the low-energy expansion, while 
“leading-order OPE” refers to the high-energy expansion. 
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out what the THS result implies for the pion distribution 
function [18, 19, 20, 60]. 


Acknowledgment 

T.H. would like to thank for the hospitality of Uppsala 
University and to all his colleagues there who made a very 
friendly atmosphere. We also acknowledge valuable dis¬ 
cussions with our colleagues K. Kampf and J. Novotny 
from Prague. We are indebted to M. Knecht for bringing 
[43] to our attention. 

The present work was supported by the Charles Uni¬ 
versity grants GAUK 700214 and SVV 260219/2015, by 
the Czech Science Foundation grant GACR 15-18080S and 
by the project “Study of Strongly Interacting Matter” 
(HadronPhysics3, Grant Agreement No. 283286) under the 
7th Framework Program of the EU. 


Appendix A. Explicit form factor formulae 


In Section 4 we have defined the invariant amplitude for 
the process 7r° —>• e + e _ using the expression P^,'P +e -. 
In this appendix we summarize its explicit form for var¬ 
ious approaches that lead to a rational function for the 
pion transition form factor. We make use of the standard 
Passarino-Veltman [53] scalar one-loop integrals Bq and 
Co- The only UV divergent function is then Bq. Its ex¬ 
plicit form will be given here as a reference point for the 
used notation: 


77r 2 B 0 (0, to 2 , to 2 ) = (27t) 4 /x 4 d J 

- IE + log 47T T log ^^ 


( 2 t T) d ll2_ m 2 + ie f 


(A.l) 


where we have introduced e = 2 — Note that differences 
of B o functions are finite. It will then be manifest from the 
form in which the results are presented, that the divergent 
parts cancel yielding a finite amplitude as desired. It is 
also convenient to introduce the following combination of 
the three-point scalar one-loop function Co and the Kallen 
triangle function A: 


z — y / / 2 2 *2 2 2 \ 

G 0 (m e ,m e ,TOi,m 2 ,TO e ,m 3 ) 

_ 1 \ <■ 2 „,2 „, 2 \ ri /™2 _„2 „ 2 „,2 „,2 _„ 2 \ 

= w\ A ( 1 ’ 2 ’ ^ C °^ 6 ’ e ’ 1 ’ 2 ’ e ’ m 3> ■ 


(A.2) 


For instance, we present two Cq functions that are used 
further in this appendix: 


Cq (to 2 , to 2 , M 2 ,0, to 2 , Ml) 

= ~ M 2 ) 2 Co(to 2 , to 2 , M 2 ,0, ml, Ml) , 


(A.3) 


Cq {ml , m e , M 2 , M Vl , m e , M Vi ) 

= -(4M 2 x - Ml) C 0 (to 2 , to 2 , Ml Ml,ml M 2 J . 

(A.4) 

Now we can finally list the results for P n o_^ e + e - as ob¬ 
tained from the various approaches to the pion transition 
form factor: 


pVMD 
■*71-°—)-e+e - 


4 

e m e 
167T 2 


Nr. 


12t r 2 F 


2 M\ 


Vi 


[B 0 (to 2 , to 2 , M^ x ) - Bq{ 0, ml ml) - 2] 


M 2 

+ Cq(to 2 , to 2 , M 2 ,0, to 2 , 0) 

- 2 Cq {ml ml Ml 0, ml Ml) 

+ Co(to 2 ,to 2 ,M 2 ,M^ i ,to 2 ,M^ i )| , 
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